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(54) METHOD AND DEVICE FOR COLLATING PICTURE 

(57)Abstract: 

PROBLEM TO BE SOLVED: To provide a method and device for collating a 
picture in which a necessary memory capacity can be reducedand complicate 
picture collation can be quickly operated at the time of collating a rotatedparallel 
movedand enlarged (reduced) picture. 

SOLUTION: A 0-(q) plane preparing part 12 prepares hi (0p) and H1 (0p) for an 
input picture by using an edge detecting part 13Hough transformation processing 
part 14Fourier transformation processing part 15and logarithmic coordinate 
transformation processing part 19and prepares H1 (0q) by performing the 
logarithmic coordinate transformation of those. Thena moving amount calculating 
part 16 compares this H1 (0q) with HO (0q) stored in a reference data storing part 
17calculates a rotational angle and an enlargement/reduction ratio in a 0-(q) 
plane leveland calculates a parallel moving amount in a 0-p plane level based on 
the pertinent rotational angle and enlargement/reduction ratio. 



CLAIMS 



[Claim(s)] 

[Claim 1]ln an image collating method which computes the angle of rotation psi of 
an inputted imagethe amount of parallel translation (xdeltaydelta)and the rate s of 
zooming to an image comparisonrespectivelyand was made to perform image 
collation of an image comparison and an inputted image based on this computed 
resultWhile performing theta-rho Hough transformation to a step which generates 
an edge direction picture about each of said image comparison and said inputted 
imageand an edge direction picture about said image comparison and generating 
reference theta-rho plane dataA step which performs theta-rho Hough 
transformation to an edge direction picture about said inputted imageand 
generates input theta-rho plane dataWhile carrying out the Fourier transform of 
the reference theta-rho plane data about said image comparison and generating 
reference theta-p plane dataEach of a step which carries out the Fourier 
transform of the input theta-rho plane data about said inputted imageand 
generates input theta-p plane datap axis of said input theta-p flat surfaceand p 
axis of said reference theta-p flat surface by carrying out logarithmic-coordinates 
conversion on q axisA step which generates the reference theta-q plane data 
concerned and input theta-q plane data so that an input theta-q flat surface may 
shift only quantity according to said rate s of zooming to q shaft orientations to a 
reference theta-q flat surfaceWhile computing a shift amount of theta shaft 
orientations of an input theta-q flat surface over a reference theta-q flat surface 
as said angle of rotation psi by comparing said reference theta-q plane data and 
input theta-q plane dataand carrying out two-dimensional image matching 
processingA step which computes a shift amount of q shaft orientations of an 
input theta-q flat surface over a reference theta-q flat surface as said rate s of 
zoomingA shift amount of theta shaft orientations of said input theta-rho flat 
surface and a rate of zooming of rho shaft orientations to said reference theta-rho 
flat surface are amended using said computed angle of rotation psi and the rate s 



of zoomingA step which computes said amount of parallel translation 
(xdeltaydelta) by calculating a cross correlation function in each theta of this 
amended reference theta-rho plane data and said input theta-rho plane dataand 
performing reverse Hough transformationand an image collating method which it 
had. 

[Claim 2]ln an image collating unit which computes the angle of rotation psi of an 
inputted imagethe amount of parallel translation (xdeltaydelta)and the rate s of 
zooming to an image comparisonrespectivelyand was made to perform image 
collation of an image comparison and an inputted image based on this computed 
resultWhile performing theta-rho Hough transformation to an edge image 
creating means which generates an edge direction picture about each of said 
image comparison and said inputted imageand an edge direction picture about 
said image comparison and generating reference theta-rho plane dataA Hough 
transforming means which performs theta-rho Hough transformation to an edge 
direction picture about said inputted imageand generates input theta-rho plane 
dataWhile carrying out the Fourier transform of the reference theta-rho plane 
data about said image comparison and generating reference theta-p plane 
dataThe Fourier conversion method which carries out the Fourier transform of 
the input theta-rho plane data about said inputted imageand generates input 
theta-p plane dataEach of p axis of said input theta-p flat surfaceand p axis of 
said reference theta-p flat surface by carrying out logarithmic-coordinates 
conversion on q axisA logarithmic-coordinates conversion method which 
generates the reference theta-q plane data concerned and input theta-q plane 
data so that an input theta-q flat surface may shift only quantity according to said 
rate s of zooming to q shaft orientations to a reference theta-q flat surfaceWhile 
computing a shift amount of theta shaft orientations of an input theta-q flat 
surface over a reference theta-q flat surface as said angle of rotation psi by 
comparing said reference theta-q plane data and input theta-q plane dataand 
carrying out two-dimensional image matching processingAn angle of rotation and 
a rate calculating means of zooming which compute a shift amount of q shaft 



orientations of an input theta-q flat surface over a reference theta-q flat surface 
as said rate s of zoomingA shift amount of theta shaft orientations of said input 
theta-rho flat surface and a rate of zooming of rho shaft orientations to said 
reference theta-rho flat surface are amended using said computed angle of 
rotation psi and the rate s of zoomingThe amount calculating means of parallel 
translation which computes said amount of parallel translation (xdeltaydelta) by 
calculating a cross correlation function in each theta of this amended reference 
theta-rho plane data and said input theta-rho plane dataand performing reverse 
Hough transformationand an image collating unit which it had. 
[Claim 3]The image collating unit according to claim 2 further provided with a 
memory measure which memorizes data about said image comparison 
beforehandand an input means which inputs said inputted image from the 
exterior. 



DETAILED DESCRIPTION 

[Detailed Description of the Invention] 
[0001] 

[Field of the lnvention]This invention about the image collating method and 
device which perform image collation of the image comparison in which 
patternssuch as a figure used as a standardwere picturizedand the inputted 
image in which the figure to be examined etc. were picturizedEspeciallyan 
inputted image rotates to an image comparisonparallel translation is carried 
outand when zooming is being carried outthe angle of rotationthe amount of 
parallel translationand the rate of zooming to an image comparison of this 
inputted image are computed at high speedand it is related with the image 
collating method and device which compare both images based on this 
computed result. 
[0002] 



[Description of the Prior Art]ln order to compare the image comparison registered 
beforehand and the inputted image inputted from the picture input device etc.it is 
necessary to compute between picturesand a relative angle of rotationthe 
amount of parallel translation and the rate of zooming. In order to compute 
thesethe technique of generalization Hough (Hough) conversion is used widely. 
[0003]There are some which were indicated by JP62-77689A as this kind of 
artand this gazette is indicating the generalization Hough transformation circuit 
which realized the circuit which can perform generalization Hough transformation 
at high speed by hardware. 

[0004]This generalization Hough transformation is an image processing 
technique for detecting arbitrary figures out of a pictureand is a method voted for 
the 4-dimensional parameter space of an angle of rotationthe amount of parallel 
translationand the rate of zooming. 

[0005]That isin generalization Hough transformationwhen performing rotation and 
collation of a picture which carried out parallel translation and carried out 
zoominga 4-dimensional memory is prepared for the parameter space 
corresponding to rotationparallel translationand zooming. And about all the 
combination between the edge points of an image comparison and an inputted 
imageit asks for the parameter of rotationparallel translationand zoomingand 
processing voted for parameter space is performed. 
[0006] 

[Problem(s) to be Solved by the lnvention]Thusaccording to generalization 
Hough transformationthat it should correspond to the 4-dimensional parameter 
space corresponding to rotationparallel translationand zoomingsince it is 
necessary to prepare a 4-dimensional memorythe problem that memory space 
becomes huge and cost will become huge invites. 

[0007]For examplein order to ask for parallel translation with resolution of 1 pixel 
and to search for an angle of rotation with resolution of 1 timewhen image size is 
256 pixels x 256 pixelsseveral 100 megabytes of memory space is required of 
parameter space. 



[0008]Since he is trying to ask for a rotationthe amount of parallel translationand 
the rate of zooming in this generalization Hough transformation from all the 
combination between each edge point of an inputted image and an image 
comparisonThe problem that the combination between edge points increases 
with the increase in the edge number of an inputted image and an image 
comparisonand processing time also increases by this invites. That isthe problem 
that collation of a complicated figure with many edge numbersetc. takes time 
occurs. 

[0009]lt is made in view of such the actual conditionrotateand this invention 
carries out parallel translationWhen you compare the picture which carried out 
zoominglet it be a solution technical problem to enable it to process promptly 
even if it is a case where the picture by which memory space could be made 
small and the complicated figure etc. were picturized is compared. 
[0010] 

[Means for Solving the Problem]Thenthe angle of rotation psi of an inputted 
image [ as opposed to an image comparison by main invention of this invention ]. 
In an image collating method which computes the amount of parallel translation 
(xdeltaydelta)and the rate s of zoomingrespectivelyand was made to perform 
image collation of an image comparison and an inputted image based on this 
computed resultWhile performing theta-rho Hough transformation to a step which 
generates an edge direction picture about each of said image comparison and 
said inputted imageand an edge direction picture about said image comparison 
and generating reference theta-rho plane dataA step which performs theta-rho 
Hough transformation to an edge direction picture about said inputted imageand 
generates input theta-rho plane dataWhile carrying out the Fourier transform of 
the reference theta-rho plane data about said image comparison and generating 
reference theta-p plane dataEach of a step which carries out the Fourier 
transform of the input theta-rho plane data about said inputted imageand 
generates input theta-p plane datap axis of said input theta-p flat surfaceand p 
axis of said reference theta-p flat surface by carrying out logarithmic-coordinates 



conversion on q axisA step which generates the reference theta-q plane data 
concerned and input theta-q plane data so that an input theta-q flat surface may 
shift only quantity according to said rate s of zooming to q shaft orientations to a 
reference theta-q flat surfaceWhile computing a shift amount of theta shaft 
orientations of an input theta-q flat surface over a reference theta-q flat surface 
as said angle of rotation psi by comparing said reference theta-q plane data and 
input theta-q plane dataand carrying out two-dimensional image matching 
processingA step which computes a shift amount of q shaft orientations of an 
input theta-q flat surface over a reference theta-q flat surface as said rate s of 
zoomingA shift amount of theta shaft orientations of said input theta-rho flat 
surface and a rate of zooming of rho shaft orientations to said reference theta-rho 
flat surface are amended using said computed angle of rotation psi and the rate s 
of zoomingln each theta of this amended reference theta-rho plane data and said 
input theta-rho plane datahe calculates a cross correlation functionand is trying to 
have with a step which computes said amount of parallel translation 
(xdeltaydelta) by performing reverse Hough transformation. 
[001 1]He is trying to have further a memory measure which memorizes data 
about said image comparison beforehandand an input means which inputs said 
inputted image from the exterior in another invention of this invention in the same 
image collating unit as an image collating method mentioned above. 
[0012] 

[Embodiment of the lnvention]Hereafteran embodiment of the invention is 
described with reference to drawings. 

[0013] Drawing 1 is a figure showing the composition of the image collating unit 
10 used by this embodiment. 

[0014]As shown in drawing 1 this image collating unit 10From the image 
comparison and inputted image which were prepared beforehandfurthercarry out 
the Fourier transform of the theta-rho plane data obtained by Hough 
transformationand Theta-p plane dataWhile computing theta-q plane data and 
asking for an angle of rotation and the rate of zooming with this theta-q plane 



levelit constitutes so that the amount of parallel translation may be calculated 
with a theta-rho plane level using the angle of rotation and the rate of zooming for 
which it asked. 

[0015]This image collating unit 10 specifically The image input part 11 as an 
input meansThe theta-q flat-surface preparing part 12 and the edge detection 
section 13 as an edge image creating meansThe Hough transformation treating 
part 14 as a Hough transforming meansand the Fourier conversion process part 
15 as a Fourier conversion methodlt consists of the movement magnitude 
calculation part 16 as an angle of rotation and the rate calculating means of 
zoomingor an amount calculating means of parallel translationthe referred data 
storage parts store 17 as a memory measurethe reverse Hough transformation 
treating part 18and the logarithmic-coordinates conversion process part 19 as a 
logarithmic-coordinates conversion method. 

[0016]Said edge image creating meansa Hough transforming meansthe Fourier 
conversion methodand a logarithmic-coordinates conversion method contain the 
theta-q flat-surface preparing part 12 respectively. 

[0017]Namelyin the theta-q flat-surface preparing part 12when creating an edge 
direction picturethe edge detection section 13 is usedWhen creating theta-rho 
plane datathe Hough transformation treating part 14 is usedwhen creating theta- 
p plane datathe Fourier conversion process part 15 is usedand also a 
logarithmic-coordinates conversion process is performed using the logarithmic- 
coordinates conversion process part 19and theta-p plane data is changed into 
theta-q plane data. 

[0018]The image input part 11 is a treating part which inputs the reference 
concentration picture R (xy) prepared beforehand and the inputting density 
picture I (xy)and outputs each inputted picture to the theta-q flat-surface 
preparing part 12. 

[0019]The edge detection section 13the Hough transformation treating part 14the 
Fourier conversion process part 15and the logarithmic-coordinates conversion 
process part 19 are used for the theta-q flat-surface preparing part 12An edge 



direction picturetheta-rho plane datatheta-p plane dataand theta-q plane data are 
created one by one from the inputted pictureand theta-rho plane data and theta-q 
plane data are outputted to the movement magnitude calculation part 16. 
[0020]An edge direction picture is specifically created from the picture inputted 
using the edge detection section 13and Hough transformation of the edge 
direction picture is carried out using the Hough transformation treating part 14 
after that. The Fourier transform of the theta-rho plane data obtained by this 
Hough transformation is carried out using the Fourier conversion process part 
15theta-p plane data is createdlogarithmic-coordinates conversion of the p axis is 
further carried out at q axis using the logarithmic-coordinates conversion process 
part 19and theta-q plane data is created. 

[0021]Namelywhen the inputted picture is the reference concentration picture R 
(xy) in this theta-q flat-surface preparing part 12. this ~ the case where the 
picture which created the reference theta-q plane data corresponding to R 
(xy)and was inputted is the inputting density picture I (xy) ~ this -- the input theta- 
q plane data corresponding to I (xy) is created. 

[0022]For examplewhen the inputted picture is R (xy)reference edge direction 
picture Re (xy) is first created using the edge detection section 13next Hough 
transformation of this Re (xy) is carried out by the Hough transformation treating 
part 14and the reference theta-rho plane data hO (thetarho) is created. The 
Fourier transform of this hO (thetarho) is carried out using the Fourier conversion 
process part 15the reference theta-p plane data HO (thetap) is 
createdfurtherlogarithmic-coordinates conversion of this HO (thetap) is carried out 
by the logarithmic-coordinates conversion process part 19and the reference 
theta-q plane data HO (thetaq) is created. 

[0023]The movement magnitude calculation part 16 which received hO (thetarho) 
and HO (thetaq) is stored in the referred data storage parts store 17 by using this 
hO (thetarho) and HO (thetaq) as referred data. 

[0024]The edge detection section 13 applies a gauss Laplacian filter to a 
concentration pictureWhile detecting an edge point in quest of the point (it is 



called the following "zero crossing point".) that the numerals of the 
positive/negative of a gauss Laplacian filter change to x shaft orientations or y 
shaft orientationsThe edge intensity Em and edge direction Etheta are calculated 
aiming at removal of a noise by applying the Sobel operator to the position of this 
zero crossing pointThis edge intensity Em is a treating part which stores the edge 
direction Etheta in an edge direction picture on condition that it is more than a 
predetermined threshold. 

[0025]Specifically by this inventionle (xy) corresponding to Re (xy) and I (xy) 
corresponding to R (xy) is created using the edge detecting method based on the 
zero crossing point indicated by JP5-151352A. 

[0026]The Hough transformation treating part 14 is a treating part which carries 
out Hough transformation of the edge direction picture which the edge detection 
section 13 createdand creates theta-rho plane dataandspecificallycreates hi 
(thetarho) corresponding to hO (thetarho) and le (xy) corresponding to Re (xy). 
[0027]The Fourier conversion process part 15 is a treating part which carries out 
the Fourier transform of the theta-rho plane data which the Hough transformation 
treating part 14 createdand creates theta-p plane dataandspecificallycreates H1 
(thetap) corresponding to HO (thetap) and hi corresponding to hO (thetarho) 
(thetarho). 

[0028]The logarithmic-coordinates conversion process part 19 is a treating part 
which carries out the logarithmic-coordinates conversion process of the p axis of 
the theta-p flat surface which the Fourier conversion process part 15 created to a 
logarithmic-coordinates axis q axisandspecificallycreates H1 (thetaq) 
corresponding to HO (thetaq) and H1 corresponding to HO (thetap) (thetap). 
[0029]The input data hi (thetarho) and H1 (thetaq) inputted from the theta-q flat- 
surface preparing part 12 and the referred data hO (thetarho) and HO (thetaq) 
stored in the referred data storage parts store 17 are used for the movement 
magnitude calculation part 16lt is a treating part which computes the angle of 
rotationthe amount of parallel translationand the rate of zooming to the figure 
(shape pattern used as a standard) included in the reference concentration 



picture R (xy) of the figure (thing to be examined) included in the inputting density 
picture I (xy)and outputs this. 

[0030]Specificallythis movement magnitude calculation part 16 computes the 
amount of parallel translation with the theta-rho plane level amended by these 
angles of rotation and the rate of zooming while computing an angle of rotation 
and the rate of zooming with the theta-q plane level which carried out the Fourier 
transform of the theta-rho flat surface which performed Hough transformation 
further. Thusthe reason for having carried out the Fourier transform of the theta- 
rho flat surface furtherand having computed the amount of parallel translation at 
the amended theta-rho flat surfaceFirstit is for asking for an angle of rotation and 
the rate of zooming promptlyafter removing the influence of parallel 
translationand calculating the amount of parallel translation promptlyafter 
removing the influence of these angles of rotation and the rate of zooming after 
that. 

[0031]When this movement magnitude calculation part 16 receives referred data 
(thetaq)i.e.hO (thetarho) and HOfrom the theta-q flat-surface preparing part 
12processing which memorizes this referred data to the referred data storage 
parts store 17 is performed. For this reasonwhen the theta-q flat-surface 
preparing part 12 outputs referred data to the movement magnitude calculation 
part 16he is trying to give the identification flag etc. which show that the data to 
output is referred data to the referred data concernedalthough explanation of 
those details is omitted. 

[0032]The referred data storage parts store 17 is a storage parts store which 
memorizes the theta-rho plane data hO (thetarho) and the theta-q plane data HO 
(thetaq) corresponding to the reference concentration picture R (xy) as referred 
dataand this referred data is accessed by the movement magnitude calculation 
part 16. 

[0033]The reverse Hough transformation treating part 18 is a treating part used 
when the movement magnitude calculation part 16 computes the amount of 
parallel translationand performs reverse Hough transformation about rho cross 



correlation picture Crho (thetarho) which specifically memorized the correlation 
coefficient which this movement magnitude calculation part 16 computed 
between the reference theta-rho flat surface and the input theta-rho flat surface. 
[0034]Thuswhen this image collating unit 10 compares the reference 
concentration picture R (xy) and the inputting density picture I (xy)While carrying 
out the Fourier transform and asking for an angle of rotation and the rate of 
zooming with a theta-q plane level after performing Hough 
transformationrespectivelybased on this angle of rotation and the rate of 
zoomingthe amount of parallel translation is computed with a theta-rho plane 
level. 

[0035]Nextthe procedure of the described image collating unit 10 is explained. 
Howeverreferred data shall already be ending with setting out here at the 
referred data storage parts store 10. 

[0036] Drawing 2 is a flow chart which shows the procedure of the image collating 
unit 10 shown in drawing 1 . 

[0037]As shown in the drawing 2if the inputting density picture I (xy) is inputtedas 
for this image collating unit 10the theta-q flat-surface preparing part 12 will create 
the edge direction picture le (xy) from I (xy) using the edge detection section 13 
(Step 201). 

[0038]And after the theta-q flat-surface preparing part's 12 carrying out Hough 
transformation of the le (xy) using the Hough transformation treating part 14 and 
creating the theta-rho plane data hi (thetarho) (Step 202)The Fourier transform 
of hi (thetarho) is further carried out using the Fourier converter 15and the theta- 
p plane data H1 (thetap) is created (Step 203). 

[0039]The theta-q flat-surface preparing part 12 carries out logarithmic- 
coordinates conversion of the p axis of a theta-p flat surface (thetap) at the 
logarithmic-coordinates axis of q axis using the logarithmic-coordinates 
conversion process part 19and creates the theta-q plane data H1 (thetaq) (Step 
204). 

[0040]And the movement magnitude calculation part 16 reads the referred data 



HO (thetaq) after the same logarithmic-coordinates conversion that was created 
beforehand and was beforehand memorized to the referred data storage parts 
store 17 from the referred data storage parts store 17 concerned (Step 205). 
[0041]Subsequentlythese two-dimensional correlation coefficients Cr (thetaq) are 
calculated using the Fourier logarithmic-coordinates resolution picture H1 
(thetaq) created at Step 204and the Fourier logarithmic-coordinates resolution 
picture HO (thetaq) read from the referred data storage parts store 17 at Step 205 
(Step 206). 

[0042]And this two-dimensional correlation coefficient Cr (thetaq) asks for the 
angle of rotation psi and the rate s of zooming from the position of theta and q 
used as the maximum (Step 207). 

[0043]Nextthis movement magnitude calculation part 16 amends the shift amount 
of theta shaft orientations of an input theta-rho flat surface and the rate of 
zooming of rho shaft orientations to a reference theta-rho flat surface based on 
these angles of rotation psi and the rate s of zoomingCalculate a normalization 
correlation coefficientshifting to rho shaft orientations about each theta of this 
amended reference theta-rho plane data hO (thetarho) and the input theta-rho 
plane data hi (thetarho)and the position of the maximum of that correlation 
coefficient is detectedrho cross correlation picture Crho (thetarho) which 
memorized the correlation coefficient is created (Step 208). 
[0044]And the movement magnitude calculation part 16 carries out reverse 
Hough transformation of rho cross correlation picture Crho (thetarho) using the 
reverse Hough transformation treating part 18creates the reverse Hough 
transformation picture Inv (xy)and computes the amount of parallel translation 
(xdeltaydelta) in quest of the position of the maximum (Step 209). 
[0045]While asking for the angle of rotation psi and the rate s of zooming with a 
theta-q plane level by performing a series of above-mentioned processingsit 
becomes possible to calculate the amount of parallel translation (xdeltaydelta) 
with a theta-rho plane level based on this angle of rotation psi and the rate s of 
zooming. In the above-mentioned procedurealthough the explanation about the 



creation procedure of Re (xy)hO (thetarho)HO (thetap)and HO (thetaq) was 
omittedlt can create by performing Steps 201-204 about these as well as the 
case of the inputting density picture I (xy). 

[0046]Nextthe creation procedure of the edge direction picture shown in Step 201 
of drawing 2 is explained concretely. 

[0047] Drawjng 3 is a flow chart which shows the creation procedure of the edge 
direction picture of Step 201 of drawing 2. 

[0048]By the edge detection section 13as shown in drawing Sin order to ask for a 
zero crossing point firstthe gauss Laplacian filter shown in a following formula to 
the inputting density picture I (xy) is appliedand Laplacian picture Ig (xy) is 
created (Step 301). 
[0049] 

And it is checked whether the noticed picture element of Ig (xy) is negativeand 
the pixel value of at least one pixel is positive among the about four pixel (Step 
302)When these conditions are satisfiedwith the application of the differentiation 
operator of Sobel (Sobel)the edge intensity Em is computed to I (xy) (Step 303). 
[0050]And if this edge intensity Em is more than a predetermined thresholdedge 
direction Etheta will be calculatedand it will store in the edge direction picture le 
(xy)and will shift to the case of less than a threshold at processing of the 
following pixel (Step 305). 

[0051]Thussince it has the character in which the thing corresponding to the 
edge point of a noise increases when the above-mentioned zero crossing point 
tends to make the value of sigma small and it is going to detect the detailed edge 
of a pictureRemoval of the zero crossing point of a noise is aimed at by applying 
the Sobel operator to the position of the zero crossing point of this noise. 
[0052]And if processing of these steps 302-305 is repeated (Step 306) and the 
processing to all the pixels is endedmaking a noticed picture element shiftthis 
edge direction image processing will be ended. 



[0053]Thusthis edge detection section 13 is creating the edge direction picture le 
(xy) using the edge detecting method indicated by JP5-151352A and the same 
technique. 

[0054]The differentiation operator of above-mentioned Sobel comprises the mask 
operator 40 of a x directionand the mask operator 41 of a y directionas shown in 
drawing 4 ln the output from the mask operator 40 of a x directionif the output 
from Mx and the mask operator 41 of a y direction is set to Mythe edge intensity 
Em and edge direction Etheta will be computed by a following formula. 
[0055] 

Although the case where the differentiation operator of Sobel was used was 
explained hereit is also possible to apply various differentiation operatorssuch as 
Roberts (Robert) and Robinson (Robinson). 

[0056]Nextthe creation procedure of the theta-rho plane data shown in Step 202 
of drawing 2 is explained concretely. 

[0057] Drawing 5 is a flow chart which shows the creation procedure of the theta- 
rho plane data of Step 202 of drawing 2 . 

[0058]As shown in drawing 5 if the edge direction picture le (xy) is inputted (Step 
501)the Hough transformation treating part 14lt checks whether this le (xy) is an 
edge point (Step 502)and in being an edge pointit sets the angle variable theta 
as theta=E theta-delta theta (Step 503). As for this deltathetaan appropriate 
value is set up based on an experiment. 
[0059]To the next 

While computing+1 addition (vote) is carried out to hi (thetarho) corresponding to 
the computed value (Step 504). 

[0060]And the angle variable theta is ************** ec | (step 505)and as long as 
this angle variable theta is within the limits of Etheta-deltatheta to E theta+delta 
thetaprocessing of the above-mentioned steps 504 and 505 is repeated (Step 



506). 

[0061 ]And after repeating this processing about all the pixels of the edge 
direction picture le (xy) (Step 507) and finishing processing of all the pixelsthis 
theta-rho Hough transformation processing is ended. 

[0062]That iswhen le (xy) is an edge pointhl (thetarho) is created by computing 
rhodisplacing the angle variable theta from Etheta-deltatheta to 
Etheta+deltathetaand adding to hi (thetarho) corresponding to the combination 
of theta and rho +one time. 

[0063]Although here explained the case where hi (thetarho) corresponding to 
the inputting density picture I (xy) was createdit can create similarly about hO 
(thetarho) corresponding to the reference concentration picture R (xy). 
[0064]Nextthe creation procedure of the Fourier conversion process shown in 
Step 203 of drawing 2 i.e.theta-p plane datais explained concretely. 
[0065]The inputting density picture I (xy) is expanded and reduced at the rate s of 
zooming to the reference concentration picture R (xy) hereit rotates with the 
angle of rotation psi furtherand parallel translation only of the amount of parallel 
translation (xdeltaydelta) is carried out further. 

[0066]Zooming on such (xy) spacerotationand parallel translation are expressed 

with the conversion shown in a lower type on space (thetarho). 

[0067] 

It corrects 
It comes out. 

[0068]Thenin order to remove the influence of such parallel translation from 
space (thetarho)M (thetarho) corresponding to the inputting density picture I 
(xy)and hO (thetarho) corresponding to the reference concentration picture R (xy) 
-- about each. The one-dimensional Fourier transform of rho shaft orientations is 
performedthe frequency p calculates the power spectral density of the field of 



p>=0 after thatand he is trying to calculate the Fourier resolution picture H1 
(thetap) and HO (thetap). 

[0069]That isas drawing 6 is a flow chart which shows the creation procedure of 
this theta-p plane data and it is shown in the drawing 6 the Fourier conversion 
process part 15 inputs first the input theta-rho plane data hi (thetarho) which the 
Hough transformation treating part 14 created (Step 601). Subsequentlyafter 
initializing the angle variable theta to zero (Step 602)FFTi.e.Fast Fourier 
Transform performs the one-dimensional Fourier transform of rho shaft 
orientations in hi (thetarho)and the power is stored in H1 (thetap). That isthe 
frequency p calculates the power spectral density of the field of p>=0and the 
Fourier resolution picture H1 (thetap) is searched for like following the (8) type. 
[0070] 

(Step 603). 

[0071 ]And processing is endedwhen it checked whether (Step 604) and ** theta 
were less than thetamax (Step 605)and processing of Steps 603 and 604 is 
repeated when it is less than thetamax and it becomes more than thetamax 
soonafter **************j ng ^ e an g| e variable theta. 

[0072]Although here explained the case where H1 (thetap) corresponding to the 
inputting density picture I (xy) was createdit can create similarly like the above- 
mentioned (7) formula about HO (thetap) corresponding to the reference 
concentration picture R (xy). 

[0073]From the above (4)(7)and (8) typesthe relation between HO (thetap) and 
H1 (thetap) is expressed like following the (9) typeand is understood that the 
influence of parallel translation is removed from comparison with (4) types. 
[0074] 

Nextthe procedure of the logarithmic-coordinates conversion process shown in 



Step 204 of drawing 2 is explained concretely. 

[0075]Herethe reason for carrying out a logarithmic-coordinates conversion 
process is as follows. 

[0076]Namelywhen it has the rate s of zooming and the inputting density picture I 
(xy) is expanding and contracting to the reference concentration picture R (xy). 
The input Fourier resolution picture H1 (thetap) is what was reduced and 
expanded to p shaft orientations to the reference Fourier resolution picture HO 
(thetap) (for examplewhen the inputting density picture I (xy) is expanded to the 
reference concentration picture R (xy)it becomes a relation shrunken by p shaft 
orientations). 

[0077]lf asked for the rate s of zooming as it isdata processing will become 
complicated and processing will take time. 

[0078]Thenthe relation which the picture expands and contracts in p shaft 
orientations is changed into the relation in which the picture is carrying out 
parallel translation to q shaft orientations by carrying out logarithmic-coordinates 
conversion of the frequency p axis at the logarithmic-coordinates axis q axis of 
frequency. That isonly the quantity according to the rate s of zooming is made to 
shift an input theta-q flat surface to q shaft orientations to a reference theta-q flat 
surface (parallel translation). 

[0079]Thusby making it the relation which is carrying out parallel translation to q 
shaft orientationsdata processing which asks for the rate s of zooming will 
become simpleand processing time will be shortened. 
[0080]As drawing 7 is a flow chart which shows the procedure of such a 
logarithmic-coordinates conversion process and it is shown in the drawing 7 the 
logarithmic-coordinates conversion process part 19 initializes q to zero while 
initializing theta to zero first (Step 701) (Step 702). 

[0081]Subsequentlyit is p corresponding to q Following the (25) type and p=c-exp 
(q) (25) 

(howeverit asks for c from constant). That isthe coordinates position on p axis 
corresponding to the coordinates position on the logarithmic-coordinates axis q is 



searched for (Step 703). 

[0082]ln this wayif the correspondence relation between p and q becomes 
clearthe Fourier resolution picture H1 (thetap) will be changed into the 
corresponding Fourier logarithmic-coordinates resolution picture H1 (thetaq) at 
following the (26) type so that it may be shown. 
[0083]H1(thetaq)=H1 (thetap) (26) 
(Step 704) 

In order to reduce the influence of the sampling error of logarithmic-coordinates 
conversionprocessing which hangs the Hanning window on q shaft orientations is 
performed to the next. 

[0084]That isas shown in following the (27) typewhat carried out the multiplication 
of the Hanning window function W (q) to H1 (thetaq) acquired at Step 704 is set 
to new H1 (thetaq). 
[0085] 

W (q) =0.5 (1+cos (piq/qmax)) 
H1(thetaq)=W (q) (27) -H1 (thetaq) 
Howeverqmax is the maximum of q (Step 705). 

[0086]Subsequently************** q (Step 706)if q is less than the maximum qmax 
(judgment YES of Step 707)will repeat the same processing (Step 703 - Step 
706) to updated qbut. If q reaches the maximum qmax soon (judgment NO of 
Step 707)it will shift to the following step 708. 

[0087]Shortly ************** thetaif theta is less than maximum thetamax 
(judgment YES of Step 708)after making q into zero again to updated thetawill 
repeat the same processing (Step 702 - Step 707)but. If theta reaches maximum 
thetamax soon (judgment NO of Step 708)this logarithmic-coordinates 
conversion process will be terminated. 

[0088]Although here explained the case where H1 (thetaq) corresponding to the 
inputting density picture I (xy) was createdit can create similarly about HO 
(thetaq) corresponding to the reference concentration picture R (xy). 
[0089]Such acquired relation between H1 (thetaq) and HO (thetaq) is expressed 



like following the (10) type. 
[0090] 

It corrects 
It comes out. 

[0091]ThusH1 (thetaq) corresponding to the inputting density picture I (xy) is q 
shaft orientations about HO (thetaq) corresponding to the reference concentration 
picture R (xy). - It can express as what shifted only psi to lambda and theta shaft 
orientations. 

[0092]Thereforeif shift amount-lambda of these q shaft orientations is 
computableit can ask for the rate s of zooming from the relation of the above- 
mentioned (12) formulaand if shift amount psi of theta shaft orientations is 
computablethe angle of rotation psi can be searched for. 
[0093]Thenshift amount of such q shaft orientations - In order to ask for lambda 
and shift amount psi of theta shaft orientations the movement magnitude 
calculation part 16Firstthe referred data HO (thetaq) after the logarithmic- 
coordinates conversion memorized beforehand is read from the referred data 
storage parts store 17 concerned to the referred data storage parts store 17 
(Step 205)These two-dimensional correlation coefficients Cr (thetaq) are 
calculated using this read Fourier logarithmic-coordinates resolution picture HO 
(thetaq) and the Fourier logarithmic-coordinates resolution picture H1 (thetaq) 
created at Step 204 (Step 206). Thenthis two-dimensional correlation coefficient 
Cr (thetaq) asks for the angle of rotation psi and the rate s of zooming from the 
position of theta and q used as the maximum (Step 207). 

[0094]As it is a flow chart which shows the procedure of calculation processing of 
data processing of such a two-dimensional correlation coefficient Cr (thetaq) and 
the angle of rotation psiand the rate s of zooming and is shown in the drawing 
Sdrawing 8 carries out the two-dimensional Fourier transform of H1 (thetaq) 



firstand calculates F1 (uv) (Step 801). Subsequentlythe power of each ingredient 
of this F1 (uv) is normalized to 1.0and it asks for F1 phi (uv) (Step 802). 
Processing with the same said of HO (thetaq) corresponding to the reference 
concentration picture R (xy) is performedand F0 phi (uv) created as a result is 
beforehand memorized by the storage parts store. Thenthis F0 phi (uv) is read 
(Step 803). 

[0095]if it explains in more detail about these the processings of a series of -- first 
-- HO (thetaq) and H1 (thetaq) -- as the two-dimensional Fourier transform of 
each is carried out and it is shown in following (13) and (14) typeFO (uv) and F1 
(uv) are calculated. 
[0096] 

Herethe expression of relations (15) of F0 (uv) and F1 (uv) is obtained using the 

above-mentioned (10) formulaand these (13) and (14) types. 

[0097] 

Thenwhen F0 (uv) is divided into power and a phaseit is expressed like following 
the (16) typeand is expressed like following the (17) type from the above- 
mentioned (15) formula about F1 (uv). 
[0098] 

As a resultthe phase constituent of F0 (uv) and F1 (uv) is expressed like 

following (18) and (19) typerespectively. 

[0099] 



Using such acquired F0 phi (uv) and F1 phi (uv)thenFlphi (uv)As shown in 



following (20)(21)and (22) typeinverse Fourier transform of the product of 
FOphi(uv) * is carried outand the correlation value (two-dimensional correlation 
coefficient) Cr (thetaq) of HO (thetaq) and H1 (thetaq) is calculated. 
[0100] 

HoweverFOphi(uv) * is a complex conjugate of FOphi (uv) (Step 804). 
[0101]Thusit turns out that the calculated correlation value Cr (thetaq) is a delta 
function (refer to the above-mentioned (22) types). And this delta function Cr 
(thetaq) takes the maximum in the coordinates position of theta=psi and q=- 
lambda. Thenthe transformation and s=k-exp (-lambda) (28) based on the above- 
mentioned (12) formula for the rate s of zooming while detecting coordinates 
position theta=psi and q=-lambda which take the maximum from Cr (thetaq) flat 
surface (Step 805) and searching for the angle of rotation psi based on this 
detection position (psi-lambda) 
It asks for (k by constant) (Step 806). 

[0102]As mentioned abovein this drawing S in order to calculate the shift amount 
of theta shaft orientations of H1 (thetaq) to HO (thetaq)and the shift amount of q 
shaft orientationsused the Fourier phase conversion method for having used the 
phase of the two-dimensional Fourier transformbut. The operation for which it 
asks by the usual matched filter or two-dimensional correlation is also possible. 
[0103]Nextthe creation procedure of rho cross correlation picture Crho (thetarho) 
shown in Step 208 of drawing 2 is explained concretely. 
[0104] Drawing 8 is a flow chart which shows the creation procedure of rho cross 
correlation picture Crho (thetarho) of Step 208 of drawing 2. 
[0105]As shown in the drawing 9 if the movement magnitude calculation part 16 
asks for the angle of rotation psi and the rate s of zooming by processing shown 
in drawing 8 The input theta-rho plane data hi (thetarho) and the reference theta- 
rho plane data hO (thetarho) are inputtedand the shift amount of theta shaft 
orientations of an input theta-rho flat surface and the rate of zooming of rho shaft 



orientations to a reference theta-rho flat surface are amended using the angle of 
rotation psi searched for and the rate s of zooming (Step 901). 
[0106]Subsequentlythe angle variable theta is initialized and (Step 902) shifted to 
zeroand -rhomax x2 is substituted for quantity deltarho (Step 903). 
[0107]And the normalization cross correlation function of hO (thetarho) and hi 
(thetarho) amended with the angle of rotation psi and the rate s of zooming is 
calculatedit memorizes in rho cross correlation picture C (thetadelta rho) (Step 
904)and deltarho is ************** e d (Step 905). 

[0108]Thenit checks whether this deltarho is less than [ rhomax x2 ] (Step 
906)and in being less than [ rhomax x2 ]it shifts to Step 904 and repeats 
processing of the above-mentioned steps 904 and 905. 

[0109]On the other handwhen the (Step 907) this angle variable theta is less than 
thetamax after **************j ng ^ e an g| e variable thetawhen deltarho becomes 
more than rhomax x2it shifts to Step 903 (Step 908). 

[0110]That isin this movement magnitude calculation part 16about each theta of 
hO (thetarho) and hi (thetarho)a one-dimensional normalization cross correlation 
function is calculated shifting to rho shaft orientationsand rho cross correlation 
picture Crho (thetarho) is created. 

[0111]lt shifts and a one-dimensional normalization cross correlation function in 

case quantity is deltarho is computed by a following formula. 

[0112] 

Nextthe creation procedure of the reverse Hough transformation picture Inv (xy) 
shown in Step 209 of drawing 2 is explained. 

[0113] Drawing 10 is a flow chart which shows the creation procedure of the 
reverse Hough transformation picture Inv (xy) of Step 209 of drawing 2 . 
[0114]As shown in the drawing 1 0 if the movement magnitude calculation part 16 
creates rho cross correlation picture Crho (thetarho) (Step 1001)after it initializes 
the angle variable theta to zeroit will set up (Step 1002)and will set -rhomax as 



rho (Step 1003). 

[0115]And this Crho (thetarho) judges whether it is larger than Cmax which is this 
maximum of Crho (thetarho) (Step 1004)and it in being bigger than CmaxWhile 
substituting this Crho (thetarho) for Cmax and updating Cmaxrho at this time is 
substituted for rhok (Step 1005). Thereforethe value of rho in case Cmax is the 
maximum is stored in this rhok. 

[01 1 6]Nextthis rho is ************** e d (step 1 006) and it checks whether rho is 
less than rhomax (Step 1007)and in being less than rhomaxit repeats processing 
of Steps 1004-1006. 

[0117]On the other handwhen rho becomes more than rhomaxreverse Hough 
transformation of Crho (thetarhok) is carried out. That isat the reverse Hough 
transformation treating part 18it is a point (thetarhok) on Crho (thetarho) 

Processing changed into the straight line shown by a ** type is performed. 
[0118]And the value of Crho (thetarhok) is added on straight-line y=-(1/tantheta) 
x+rho k/sintheta on the Inv (xy) flat surface after this reverse Hough 
transformation (Step 1008). 

[01 1 9]And after ************** ing the ang | e variable thetait checks whether (Step 
1009) and ** theta are less than thetamax (Step 1010)and in being less than 
thetamaxit shifts to Step 1003. 

[0120]That isin this movement magnitude calculation part 16the position 

(thetarhok) which takes the maximum by each theta is detected from Crho 

(thetarho)reverse Hough transformation is performed in that positionand the 

reverse Hough transformation picture Inv (xy) is created. 

[0121]And the position (XmaxYmax) of the maximum of this reverse Hough 

transformation picture Inv (xy) serves as the amount of parallel translation 

(xdeltaydelta). 

[0122]Nextthe processing result at the time of applying the image collating unit 
10 concerning this invention to collation of a character is explained. 
[0123] Drawing 11 is a photograph which shows the halftone image displayed on 



a displaywhen the image collating unit 10 shown in drawing 1 is applied to 
collation of a character. 

[0124]Heredrawing 1 1 (a) is the reference concentration picture R (xy) of a 
character ("10000")and comes to show drawing 1 1 (b) reference edge direction 
picture Re (xy) which extracted the edge direction from this reference 
concentration picture R (xy). 

[0125]And if it becomes the reference theta-rho plane data hO (thetarho) which 
has a band-like pattern as shown in drawing 11 (e) when Hough transformation is 
performed to this Re (xy) and the Fourier transform of this hO (thetarho) is carried 
out furtherThe reference theta-p plane data HO (thetap) as shown in drawing 1 1 
(g) is obtained. If logarithmic-coordinates conversion of this p axis of HO (thetap) 
is furthermore carried out at q axisthe reference theta-q plane data HO (thetaq) 
as shown in drawing 1 1 (i) will be obtained. 

[0126]On the other hand drawing 1 1 (c) is the inputting density picture I (xy) of a 
character ("10000"). 

[0127]For exampleif the case where the character of the label which stick and is 
in a box is read optically is assumed in order to classify a boxas shown in this 
drawing 1 1 (c)to the character on the reference concentration picture R (xy)it 
rotatedand the read character might carry out parallel translation and may be 
expanded (reduction). 

[0128]The input edge direction picture le (xy) which extracted the edge direction 
from this inputting density picture I (xy) comes to be shown in drawing 11 (d). 
[0129]And if it becomes the input theta-rho plane data hi (thetarho) when Hough 
transformation is performed to this le (xy)as shown in drawing 1 1 (f) and the 
Fourier transform of this hi (thetarho) is carried out furtherthe input theta-p plane 
data H1 (thetap) as shown in drawing 11 (h) will be obtained. If logarithmic- 
coordinates conversion of this p axis of H1 (thetap) is furthermore carried out at q 
axisthe input theta-q plane data H1 (thetaq) as shown in drawing 1 1 (j) will be 
obtained. 

[0130]Herethe reason the wave of a belt is seen for this hi (thetarho) as 



compared with hO (thetarho) is because the attitude angle and position of a 
character differ from each other. 

[0131]Since the inputting density picture I (xy) is expanded to the reference 
concentration picture R (xy)the input Fourier resolution picture H1 (thetap) is 
understood that it has shrunken in p shaft orientations to the reference Fourier 
resolution picture HO (thetap). 

[0132]When the relation whose picture has shrunken with such p shaft 
orientations sees H1 (thetaq) and HO by which logarithmic-coordinates 
conversion was carried out (thetaq)it is understood are changed into the relation 
in which the picture is carrying out parallel translation to q shaft orientations. That 
isit turns out that that to which the picture HO (thetaq) was shifted caudad is the 
picture H1 (thetaq). 

[0133]Now drawing 11 (k) is a correlation value Cr (thetaq) flat surfaceand it turns 
out that it has a coordinates position of the maximum luminance which this 
correlation value Cr (thetaq) flat surface can distinguish from other points clearly 
(the coordinates positions which take maximum luminance are theta=psi and q=- 
lambda). 

[0134]About each theta of hi (thetarho) and hO (thetarho) drawing 1 1 (I) is a 
shown figure rho cross correlation picture Crho (thetarho) calculated while 
shifting to rho shaft orientationsand drawing 11 (m)lt is the reverse Hough 
transformation picture Inv (xy) which detected the position (thetarhok) which has 
the maximum by each theta from Crho (thetarho)and performed reverse Hough 
transformation in the position. 

[0135]As mentioned abovein this embodimentthe theta-q flat-surface preparing 
part 12 creates hi (thetarho) and H1 corresponding to an inputted image (thetap) 
using the edge detection section 13the Hough transformation treating part 14the 
Fourier conversion process part 15and the logarithmic-coordinates conversion 
process part 19While creating H1 (thetaq) which furthermore carried out 
logarithmic-coordinates conversion of thiscomparing HO (thetaq) which the 
movement magnitude calculation part 16 memorized to this H1 (thetaq) and 



referred data storage parts store 17 and asking for the angle of rotation psi and 
the rate s of zooming with a theta-q plane levelControlling memory size required 
for collationsince it constituted so that the amount of parallel translation 
(xdeltaydelta) might be computed with a theta-rho plane level based on this angle 
of rotation psi and the rate s of zoomingeven if it is a case where the edge 
number of a picture increasesthe processing time required can be reduced. 
[0136] 

[Effect of the lnvention]As explained abovewhile computing an angle of rotation 
and the rate of zooming with the theta-q plane level which carried out the Fourier 
transform of the theta-rho flat surface which performed Hough transformation 
further in this inventionSince the amount of parallel translation was computed 
with the theta-rho plane level amended by these angles of rotation and the rate of 
zoomingWhile being able to ask for an angle of rotation and the rate of zooming 
promptly after removing the influence of parallel translationafter removing the 
influence of these angles of rotation and the rate of zoomingthe amount of 
parallel translation can be calculated promptly. And a note can be made 
simultaneously and capacity can also be made small. 

[0137]For this reasonfrom generalization Hough transformationas compared with 
the case where it asks for an angle of rotationthe amount of parallel 
translationand the rate of zoomingmemory space is reduced by leaps and 
boundsand even if it is a case where the edge number of a picture 
increasesprocessing time is shortened substantially. 

[0138]When asking for an angle of rotation and the rate of zoominga frequency p 
axis by carrying out logarithmic-coordinates conversion on the logarithmic- 
coordinates axis q axis of frequency especiallySince he is trying for an input 
theta-q flat surface to shift only the quantity according to the rate s of zooming to 
q shaft orientations to a reference theta-q flat surfacedata processing which asks 
for the rate of zooming will become simpleand processing time can be shortened 
by leaps and bounds. 

[0139]The memory measure which memorizes the data about an image 



comparison beforehand in another invention of this inventionSince it was made to 
make it memorize beforehand before it had the input means which inputs an 
inputted image from the exterior and the inputted image inputted the data about 
an image comparison from the outsidelf an inputted image is inputtedit will 
become possible to compare an image comparison and an inputted image 
immediatelyand processing time can be further shortened now. 



DESCRIPTION OF DRAWINGS 

[Brief Description of the Drawings] 

[Drawing 1]Drawing 1 is a figure showing the embodiment of the image collating 
unit concerning this invention. 

[Drawing 2]Prawing 2 is a figure showing the embodiment of the image collating 

method concerning this inventionand is a flow chart which shows the procedure 

performed with the image collating unit shown in drawing 1 . 

[Drawing SjDrawing 3 is a flow chart which shows the creation procedure of the 

edge direction picture which the edge detection section shown in drawing 1 

performs. 

[Drawing 4]Drawing 4 is a figure showing the differentiation operator of Sobel. 
[Drawing 5]Drawing 5 is a flow chart which shows the creation procedure of the 
theta-rho plane data which the Hough transformation treating part shown in 
drawing 1 performs. 

[Drawing SjDrawing 8 is a flow chart which shows the creation procedure of the 
theta-p plane data which the Fourier conversion process part shown in drawing 1 
performs. 

[Drawing 7jDrawing 7 is a flow chart which shows the procedure of the 
logarithmic-coordinates conversion process which the logarithmic-coordinates 
conversion process part shown in drawing 1 performs. 

[Drawing SjDrawing 8 is a flow chart which shows the calculation procedure of an 



angle of rotation and the rate of zooming which the movement magnitude 
calculation part shown in drawing 1 follows. 

[Drawing 9]Drawing 9 is a flow chart which shows the creation procedure of rho 
cross correlation picture Crho (thetarho) which the movement magnitude 
calculation part shown in drawing 1 performs. 

[Drawing 10]Drawing 10 is a flow chart which shows the creation procedure of 

the reverse Hough transformation picture Inv (xy) which the reverse Hough 

transformation treating part shown in drawing 1 performs. 

[Drawing HjDrawing 11 (a) - drawing 11 (m) are each photograph which shows 

the halftone image displayed on a displaywhen the image collating unit shown in 

drawing 1 is applied to collation of a character. 

[Description of Notations] 

10 -- An image collating unit and 11 --An image input part and 12 -- Theta-f flat- 
surface preparing parti 3 [ -- A movement magnitude calculation parti 7 / -- A 
referred data storage parts store and 18 / ~ A reverse Hough transformation 
treating parti 9 / -- A logarithmic-coordinates conversion process part4041 / ~ 
Sobel operator ] - An edge detection section and 14 - A Hough transformation 
treating parti 5 -- The Fourier conversion process part and 16 
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tuE#P,s 0 - p ¥®lc*f T^fuEAT: 0 - p ¥®<D 0 Hi 
TJft<D->7 h«33<fclf pHi£lR]®ffi* • f§/J«ffiIE 
L, COSiE*ti/-c#fiS 0 - p¥It-^ ituEAT: 0 
- P^St 1 — 5»©*0 Jcfi^T«BBIH«t**»s a* 
/\7^^ToTtuE¥ : fT^K)M (x A, y A) ^Jtffi 

-r«¥^j»»«3tai#at^A/-cB«fis^K«o 

[«*^3] tuE#PSB#lcHr^x-^^46 
E'ft LT33 < E«#atx luEA7]B«^SP<fe U AT: 

Bo 

[0 0 0 1] 

®/^-><Wi«*ti/-c#i!SB«t«a*fftOia^^ 
Ji«**i/-c A^BftOBftfiS-g-*^ a B«fiS-&*S6Stf 
SBlcl L, #lcAT:B^#PSB#lc» LTHIeL, 
¥ff»tL ffiA-H^LTl^^lc, KA7]B«0 
#fiSB«fc*f-r«0Sftx ¥^j»B«fi<fetfffi* • f§'> 
*^iiSlc»[iJ LTs C03tail£micS^^TBB«€' 

as-g-r *b«sb-&£5£ sifsmc ir 

[0 0 0 2] 

^46SS*ti/-c#fiSB«t. B#AT:g 
B^^ 6 AT: * fife A*iB«£SB-&-f « fc 46 lc tt, B# 

-SSfb/\7 (Hough) $Sj£^-5#5£#l£<fl3^S*i 

[0 0 0 3] COfflOSffit LTm7L\mfflV& 6 2-7 
7 6 8 Q^afBlC^/Ts^tl/ctcD^Us CO^att. 

THH L/c-)Efb/\7^lHlK^F^ LTl^o 

[0004] co-ttfb/\7SSittt. fi*oia»*H 
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[0005] r^fe-gs -mti\7Wikz«&s mi, ¥ 

[0 0 0 6] 

7mmz&tu& m^ ^mm, as* • «i>Nc*usf 

[0 0 0 7] mX&s W-mtWk 1 BfR<D5Tfl?f£T\ 0 
5 6Ifx2 5 6IiiL/ct^ /O^-^^Hffllc 
[0 0 0 8] $/c, 6^-fS{b/\7^T1^ A73B 

A7:Bfs<!:#PSBfKDx -y S^ROJfiP leff o Ti 
•y s Jjm<mtt&M± L^tilCcfco TSQIIBIH & ±1* 

&0J^®^fcttBvTlH#jb\fr* t rm&K&T 

[0 0 0 9] LfcH^lCig^T^^tlfc 

[0 0 10] 

[KH^ft?)*r^/c46(D#IS] *3&E®£fc* 
38ET?tt. #SBH«fcfcf-f*A73ii«®lH]£ft^ ¥fi 
iftl (x A, y A) JJcfctfffiA • ffi/J\* s 

»aj u K3tansmic»^T#fisii«t A73B«oii 

«BB^£f55 £ 5 lc LfciMMB^KS fcfil^T. tuE# 
BP^£otT/fuEA7:Bti©^*l^*Uco^TX >y 915 
|p]B«££fiBf3X^ v 7<t, fuE#BBB«fcoi/»T® 
x >y v^ftBfilcft LT 0 - p /\7tS?fro T#BB 0 
- pW-MT-ttr^mtZtt^. fuEA73B«fco 
l^TOxy$;SiRiB«fc*fLT0-p/\7$Si*fToT 
A7] 0 - pf ET-^f T^tv MEWSR 
B«lco^TO#fiS 0 - pfffir- ££7- U XS«L 
T#!I0-p¥1t- Si ttmt* tt fcfc, i$EA7J 
Biilco^TcDA7:0-p¥S7 :, -^^7-Ux^L 
TA*e-p¥l7-**SfiEmf'^is fuEA 



7J0-p¥ScDplii SuE#18 0 - p ¥®G> p ft® 

en*, qHiic*f8tjs«$Si-r«ctic<fey. #pse- 

q¥®fc*fLTA73 0-q¥®#. qlATf^CfuEffi* 

0 - q flf- StfcJctfATJ 0 - q ¥Sx-^^J5gr 
5XT7 7"i, SuE#BS 0 - q ¥Ex-£ 4: AT: 0 - q 
flf- Sit^#-&to-&T2^7uBfi!T'y^y^jKiS 
f 3 C £ IC cfc U s #PS 0 - q ¥ElC*fr 3 AT: 0 - q ¥ 
ScD0|i*rS]CD->7 h«£fuElE]Kft0i: LTHaVTS 
iitlC, #BS0-q^MC»r^A7:0-q¥®CDq 
W5 fa(Dis7 h«£fuEffi* • ffi/J\* s t LTStii <T3 

* s £ffl^TfuK#BHe - p¥®fc*ff *fuKA*l 0 - 
p^FScD0|i7^rS]CD->7 h«fi*tf pHlSlRlOffi* • f§ 

I5IBA7: 0 - p¥Sx-^cd#0 fcte^z^mmim 

£#4k iS/\7£Sj£ftoTfuE¥ffi»»l« (xA,y 
[0 0 1 1] $/c, *£E®BU®3&ET?tt. ±5EL/cB 

vmsum t mm&m9tm£&m.fc tsi^z, iuk#sbb 
A*B«*^sp<feyA7]-r«A7]#at*sic^As<fe 

[0 0 12] 

[0 0 1 3] 0 Hi, *Hfifi®J&HT?fl3l^B«SB-SS 
B1 OODi^^tH^SS,, 

[0 0 14] 0i(c/Tvt<t:3lc, coa«fiS-&S«i o 
tt, ^46*«L/-c#fiSB«tA7]B«A^5. M7£&lc 

e fcyff6ti«0-p¥ST f -^^*eic7-yx^L 

SlC0-q¥Sx-^^»tiJU 
fe^iitlc, 5R46/c0IeS, ffiA-l§/J«ffl^T0 

[0 0 15] ^ftWlctt. COB«fiS-&S«1 Ott, A 
TJ^ISi LT©B#A7:gP 1 U> 0-q¥®f^fiESP1 
Its X -y v^B^^ISi LT0X «y v^ffigp 1 3 
/\7$Si#atLTO/\7$SiftlSSP1 At, 7- 

u x$si#at LT07- y x^sqiisp i 5 1 , 01s 
ftfi<fetfffi* • iB/j\*3tai#as^Ltt¥^j»Bi«3tai 

#at LT(D^K)M»[iJgP 1 6 <t, E«#at LTO# 
^t— $tBffff1 It, 5S!/\7$SiftlSSP1 
JS«$Si#at LT®tt»J£«$SjftyiSB 1 9tfrtt 

[0 0 16] fpj> luExy5?B«4fiE#at. /\7$Si 

#a<t, 7-yxs«#at. smsi^fi^att, # 
[0017] -rato-Sx 0-q^s^aui 2m x 
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y v^Bf^fFJST *BSfcx v SXttffltf 1 3 

[0 0 18] B«A*SP1 1 14, ^4&*fi§L/-c#SB*6 
H«R (x, y) <tA7J;lJtBfH (x, y) *A*f 

^sqiispt-s y > a* $ti&£M& e-q ^mftmm 
1 uautiTZo 

[0 0 19] 0-q¥®f^fiESP1 214, Xy$*WH»1 

tWStJSaSSiftlSSP 1 9£ffl^T, A^tl/cBfifr 
0-pfi7-^ 0-p¥S7- 
4»fi<fetf 0 - q¥Sx-^^IH^J5gU 0 -p¥l7 

[0 0 2 0] mttWICtt, Xy$*WH»1 3£ffl^TA 
* **l/-cii«fr 617 $>£|R]ii«£fMB U *®&/ \ 7 

-^*7-yx$siiQsgpi 5£fl3ivr7-yx$gj 

U 0-p¥Ex-^^JSU *5lc*f»JSa$SiJQ 
SSP1 9£ffl^Tpii£qiilCftiMI£iflLT, 0- 
q¥l7-^MtS 0 

[0 0 2 1] f&fo-Sx C®e-q¥®fMBSB1 2 7? 
14. A7J*ftfcH«##BB»6ii«R (x, y) 

RR (x, y) lc»jSr^#PS0-q^Sx 

-*£f&SU *fcA**tifc«W<A*»Sli« I 
(x, y) T?S*J§^fctt. IS I (x, y) fcfcUSf* 

A^0-q¥Sx-^^J5gr^o 
[0 0 2 2] #iJa«\ A7:£*l/cBfi#R (x, y) V 

zzmsiaz, $r««3(cx>y>>^aspi 3*«^t# 

BBxv^fiHfcRe (x, y) «ft|«U #lc, Z(D 
Re (x, y) */\7$SISQSSP1 4lcJ:y/\7$Si 
L, #PS0- p ^Sx-^hO (0, p) £fMBfS 0 
£6lc, £<DhO (0, p) *7-yx$«iiQSgpi 5 
*Jllvr7-yx**U #PS0-p¥S7 :, -^HO 
(0, p) «ft|«U *6fc, £<DH0 (0, p) £fcf 

gfcga^swiassB 1 9 ic cfc y *fa&a£« lt, mm e 

-q¥Ir- £H0 (0, q) £fMB*-*o 
[0 0 2 3] hO (0, p) SlfHO (0, q) ^gttlX 
o/cSK)M»[iJgP1 614, £<DhO (0, p) SlfHO 

(0, q) *m%T-$tLT, m^T-timm-^ i 

[0 0 2 4] xy^«aiSP1 314, JljtBfilcft LTtj 
■7X57 : 7~>7'y7-f/l/*-:&>IfflLT, 130X^7^ 

zs7>74iu$ -(om^oMm? x iteiRixtt y »£ir] 



/-r X©fifc££EI y OOI7 >>'3SJg E mStf x >y 5>#|q| 
E 0 L, ii >y $>3ffiS E m&mfe<0 L#lMiIJ}(± 
Z«tZ>Z>Z.t*:$&£ LT*©X>y$>75lR]E 0£x>y v% 

iRiii«ic«tt-r«ftissP7S«o 

[0 0 2 5] jy*TO4, *5ITO4, #l¥5-1 5 
1 3 5 2^a$glC^Tx?tl^€P^PX*lcS^<X-y 
5*l*ffl*5S*Jll»vr* R (x, y) fcfcfjSfSRe 

(x, y) SSefctf I (x, y) [<Ml&t%> I e (x, 
y) £fMBf* 0 

[0 0 2 6] /\7$SiftlSgP1 414, xy$*WH»1 3 
jJWMB L/cX y =J15 $mm:l \7%m LT 0 - p ¥®T 
-*£fMB-f*ft!3SBT?£y. jy*TO4, Re (x, 
y) |CJtJST*hO (0, p) S3<fclf I e (x, y) lc 
fcfJSf*h1 (0, p) £fMBf* 0 

[0027] 7-ui*«woawi 5i4, j\7mmmw 

SP1 4jWfiBL/-ce-p¥®7 ? -*£7-yx£JfcLT 
0 - p¥Ex-££ftjST«IBPT-£y, Jlttftlc 
14, hO (0, p) KfcfjSfSHO (0, p) e«felfh1 

(0, p) khim-shi (0, P ) zftmtZo 
[0028] tmrnumtsmimi 914, 7-yxg» 

SQIISP 1 5 jJWMB L/c 0 - p ¥®<D p ftifcjSaii 

q m^mmm^mmmt ««isspt?s y , jy*wic 

14, HO (0, p) KfcfjSfSHO (0, q) S3<fctfH1 
(0, p) fcfcfj£f*H1 (0, q) £fMBf*o 
[0 0 2 9] ^K)M»tiJgP1 614, 0-q¥®f^fiESP1 

2^6A7:?tl^A7l7 :, -^h1 (0, p) StfHI 
(0, q) t#«7-' £IE«SP1 7(C*SI«*tlfc#ftgT 

-^hO (0, p) SlfHO (0, q) ^ffll^T, A7J;1 

JSB^I (x, y) iC^ZtlZWB {^MMlR) O, 

#ps;ij*B#R (x, y) ic^$ti^0ff$ mmtr+z> 
[0030] m^mat, co»Ki«3taispi 6t4, n 

7^Sl^ofc 0 - p ¥S^? 6 lc 7- y x$«i L/c 0 

- q¥®U'W?iH]£fta3<fcifffi* • f§'j«iitij-r3 

?ftfc0-p¥EU^/UWT^ftM£StiJ-f3o 
cfcaic0-p¥S^?6lC7-yx^L, SIEL/C0 
-p¥ST-¥^K)M^»[iJL/cIIS{4, £T> ¥?t^ 
K)CDff^^^£L/c±'c:-0lES33cfcy : ffiA • ffi/J\**ffl, 

[0 0 3 1] $fc, ^(D^K)M»[iJgP1 66\ #8^^- 
*-f&:b-5h0 (0, p) SlfHO (0, q) ^0-q¥ 

/cd6, ; e<D|fflSM0Sl4*SSr^/)\ 0-q¥®f^fiESP 

1 2 6w»iM*ajgp 1 6 k*j ^TrnkT-ttfth-** 
i4, mar* x- ^ - * « c <t 
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[0 0 3 2] #PfiT~*KttSn 714, #B3££B&R 
(x, y) KfcfJS?-S0-p¥®7 ? -*hO (0, p) 
Stfe-qTSx-^HO (0, q) ^tiT-^iL 
TEBf *EttS5T?S U s frfrJMWHx-* 1*^Km» 
mgp i 6ic<fey7 , ^-tx*ti« 0 

[0 0 3 3] as/\7^SQHSP1 814, ^K)M»[iJgP1 

^ttWlctt. ^(D^K)M»[iJgPi 6##P,S0-p¥E£ 
A7J 0 - p¥®tOIHT?3tai L/cfgfU^iMB'fl L/cp 

mmmm^cp (0, P ) ^^zmj\ymm^mn 

[0 0 3 4] Z<0&ol^ COB«BS-&gB1 014, # 
S3B&B&R (x, y) <kA7J»6B«l (x, y) £ 

ss-s-f *isfc\ %ti?tu\7mtSi*fi-otcmc7- y x 

£&LT 0 - q ¥® U^/l/^iaSft fi<fetfffi* • ffi/J\* 
£;El<#>3 <!:<!: tic, RI3£fta3<fclfffi* • ffi/J\*lcS-3 
l^T 0 - p¥® U^/l/7¥ffi»KiB£J*aj-f St® 1 ?* 

[0 0 3 5] ±EB«fiS^g«1 0(DSQII#IHIC 

fMBx-^KBSB 1 0 lcR£SI*T?S« j £»Ofr«o 

[0 0 3 6] 0 2 14, 01 ICifxT BfsP,s-&gB 1 o <D$Q 
IIW&ifxT 7P-7-f- hT-S^o 
[0 0 3 7] ^0 2 ICjfxT <fc a lc, Z. OB«BS-&gB 1 
014, A7J»6B«I (x, y) ^A73**l«tx 0- 
q ¥®fMBSB 1 2 #l 7 S^ffigp 1 3£-ffll/>T I ( x , 
y) ^5iy^lRli«l e (x, y) £fMBf« (A 

7772 0 1) 0 

[0 0 3 8] ^LT, 0-q¥®fMBSB1 2tt, M7£ 
SjiQggBI 4=&ffll^T I e (x, y) £/\7£SILT0 
-p¥Ex-*h1 (0, p) *f&SLfc« (Xt77 
2 0 2) , 7-yx*«Wl 5*ffl^Th1 (0, p) 
**6tC7-yx^SiL0-p¥ST f — S» HI (0, 
p) £fMBf« 0 3) 0 

[0 0 3 9] £SlC, 0-q¥®fMBSB1 214, ftife* 
WSSiftlSSPI 9^ffl^T0-p¥S (0, p) ©pli 

qHit^-5*fajs«iiiic*fajs«$SiLT. 0-q 

¥Et-*H1 (0, q) £fMBf* U7772 0 

4) . 

[0 0 4 0] -f-LT, ^ftmittiJBP 1 614, WFJ^L 
#Pf 7-*EHSB 1 7 ic^&E'fl LTte^tcmW&im 

mmmmfeommT-tHo (0, q ) ^mmm?- 

*E«SB1 7frSfI^tiJT (X77 72 0 5) „ 
[0 0 4 1] X7772 0 4Tft)|Jtlft7- 

Uxfcf»j£«£SjB«H1 (0, q) 77772 0 
5 T?*^— *K«W 1 7 #SR*ffi*Jfcfc7- y IM 

mmmmmmmo (0, q) <t£ffli^T, <me>(D2 

&7cfflB8ffi»Cr (0, q) fctWTS (X77 72 0 



6) 0 

[0 0 4 2] ^LT, C02:&7cffll!Bffi»Cr (0, 
q) 0, qOffi«<feyiHlSft^ ffi*-f§ 

'J^»s*«a6* (X77 72 0 7) o 

[0 0 4 3] Mts ^(D^K)M»[iJgP1 614, Cti50 

n> ts&vt&x • m^m s k*^t#bs 0 - p ¥® 

lc*fT3A7: 0 - p¥S00li75^cD->7 h«fi<fetf p 
ii75ft<Dffi* • f§'J««IE L, £<D«IE L/c#fS 0 - 
p¥Sr- St hO (0, p) <tA7:0-p¥E7-£h1 

(0, p) (D : &0lco^Tpli75[R]lc->7 hL&#SIE 
&fbil!Bffi»£W-J* U *Offll!Bffi»®»*fiI®ffiB«- 
IttiJLT, fflBB«a*E«L/-cpiiSiiBBB«Cp 

(0, p) £fMBf* (7,77 72 0 8) 0 

[0 0 4 4] -f-LT, ^ftmSttiSP 1 614, i8!/\7£Sj 
ftlSSPI 8*ffl^TpfflSiBBB«Cp (0, p) «S» 
/\7^LT3S/\7^B# I nv (x, y) Sfftg 
U *®B*fiI®ffiB£$4&T¥ffi&KlB (xA,y 
A) =&lttiJT3 (777 72 0 9) 0 

[0 0 4 5] ±E-iiOftlS^5CtlC<feyx 0-q 

-p¥il/^R-¥fTili (x A, y A) ^46^^ 
t^RTlltSSo ±ESQH#IHlcJ5^Tl4s Re 

(x, y) , hO (0, p) , HO (0, p) , HO 

(0, q ) (D\mmm^^T(Dmw*^Ltci5\ z. 

*i6tco^T*>, A7J»6E«l (x. y) o>m;m 
«IC77772 0 1-2 0 4£5Sftf*C£:fc:<fcyftoB 

[0 0 4 6] 02(7)77772 0 1 KjjVTxy-7 

SlRlB«0^fiE#lllco^T^ft«llcWWr«o 
[0 0 4 7] 0 314, El 207x7 72 0 1 0X7-775 
[R]B#0^J5K#|iI : &/Txr7P-7l'- h^SSo 

[0 0 4 8] 03 ictxT cfc a ic, x 7 s;«ajgp 1 3 t- 

14, $rftWlc€P7P7*^46^/c46lc, A7:;lJS 
BfH (x, y) lc»LT^lOTs^777 77->7 7 
V7Y/P'S'— £jffflU 77x>7>If I g (x, 
y) ZftmtZ (777 73 01), 
[0 0 4 9] 

V 2 GC,,,) = ^(l-^!) e -^ (1) 

■f-LTs I g (x, y) oaSB*^^ ^© 
43£fi|0BJR<Od'5^ <<t«6» 1 O(DB*0B*fi^IE 
T»»S6^«iiRL (7x7 73 0 2) % 

ti>m±L?zm-£im, \ (x, y ) icwltv-^/p 

(Sobel ) (D^7|-^U-^^3ifflLTX7v >, ?tJSEm 
^»tijr^ (77773 0 3) o 

[0 0 5 0] ^-LTs C(DX7v >, ?tJSEm^mS®L* 
l>filJJ(±T-«*ltf> X 7 -775ft E 0 *W-Jt LTX 7 -775 
ftB^I e (x, y) icftttU L*^fifi*3l©»&lc 
14^ *B*CDSQIIlC^7f§ (777 73 0 5) 0 
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[0 0 5 1] C0<fc5lz\ ±K-tfn*PXj5tt. o(Om 

*/h*< LTmmvwmKxv'jz&iiiL&otTzt 

[0 0 5 2] -f-LTs aSB**i»^**S^5A^« 
Xr^3 0 2~3 0 5 OiQg^SK y >g L Uf7^3 
0 6) , ^TcDH«lc»rWI^H7L/c^6«; C 

[0 0 5 3] ^(DcfcaiC, COXy^ffigpi 3tt, ft 
HT5-1 5 1 3 5 2^a$glcr^Tx?tl/cX-y^tiJ* 
»4:BMW>#a*ffl^Txy$;SiRili«l e (x, y) 

[0 0 5 4] ±IBV-^/KD^^U-^^ 
04tC^-TJ:3^ x£|r]0TX-7:>|-^U~$14 0<ty 
1ofa(D^X<7*«.ls-$4 1 frS^J^tU x£|r]®T 
X^^U— £4Ofr60iB2)£Mx, y£|R]®TX? 
*^U-^4 1 A^50ai**Myfr«tx X-yv^tj* 

[0 0 5 5] 

E B = tan" 1 ^) (2) 

•&lC"O^TI/^L/c6\ P/t-y (Robert) ^>Ptf>V 
V (Robinson) tW)«Wt»*^U-< ££i!fflr3C£: 

[0 0 5 6] m2(DXTv72 0 2fcmt9-p 

[0 0 5 7] H5ttx E12cDXx i y^2 0 2(D0-p¥ 
If- £ Qftti&mtemty P + - h T?S 

[0 0 5 8] 0 5lC,Tx1-cfcaiC, /\7$SlftlSSP1 4 
tt, iyi^Ml e (x, y) #A7J?*l3<!: (X 
7775 0 1) , Z(D\e (x, y) tfxy 
frS#£5IKL Uf775 0 2) s xy 

0 = E 0-A0 

lcfS£T3 Uf775 0 3) 0 ZCDA91Z. H 

[0 0 5 9] :Mc, 

P = xcasO + y&n& 
(Ee-Ad<0<E e +Ae) (3) 

^saj-f^itic, *©3taifiiic*fjs-r« h i (0, 
p) tc+ijjp» (fss) r« (xf775 0 4) 0 
[0 0 6 0] -f-LTs ftj^»e*-r>^y^>hu 

(7.7^75 0 5) % RftS£R0#E0-A0#SE 
0 + A0©l5HrtT»»Siayx ±IBX7 1 'y75 0 4&tf 



5 0 5©jflS*l»yig-r (Xx77"5 0 6) 0 
[0 0 6 1] ^-LTs frfr«I=&x-y v>£ftlHi I e 
(x, y) ©£HJgfcOl^T«iyjgL Uj775 0 

7) , ^U^CDSQII^A/cBi*^ C00-p/\7£ 

[0 0 6 2] T&totS, I e (x, y) #Xy5W* 
ftJS$»0*E0-A0A^5E0 + A0$ 
TtfuftOOp^IffiU *O0StfpOffl*-&*5* 
KfcfJSfShl (0, p) lc+iiP3tLT^<Ci:tcJ: 
>A hi (0, p) *f&*LTl^6o 

[0 0 6 3] CCT?ttA*l»6ii«l (x, y) 

fcfcfj£f*h1 (0, p) £fMBf*»&fc:o^TKE 
L/c6\ #P,S;lJtiHiR ( X , y) [ZtjfctZhO 

(0, p) lco^Tt»Mt§Ci^t§o 

[0 0 6 4] 0207.^72 0 3lCifxT7-y 

xSSIiQS. OS y 0 - p¥®T-*®fMB#|lIfc:o^ 

[0 0 6 5] A7J;1JSH#I (x, y) » s # 

BftX&Sft-R (x, y) fcfcfLT. ffi* • ffi/h* s T?Jfi 

(x A, y A) Ett¥fj|«bStlT^*o 
[0 0 6 6] EdLfe (x,y) SIH±T?Offi* • « 
'K (0,p) ffilH±7ttTaKw 

[0 0 6 7] 

mm - ^(^-y /" fCQ f" K) ) w 

fc/cU 

* - V^nd (5) 

K - tan " 1( ^! < 6 > 

[0 0 6 8] (0,p) ffiMfrS, £?Lfc¥ 

fit$W}(D§m%:Vk£tZ>tclsf)fc, A7J;1JSH# I (x, 
y) KfcfJSfShl (0, p) e*tf#fiS»Sli«R 

(x, y) KfcfjSfShO (0, p) WftUcoi^ 

pa*iRioi*7c07-yx$si^i\ ^cDfgji 

S» p # p ^ 0 (DfIi^cD/\°7-A^<7 h/l/»6£trJ* L 
T7-yx$«ili«H1 (0,p) > HO (0, p) *m 

[0 0 6 9] T^*)"^ me-pTffir 
-^(DftjSW&ifx-f 7P-t>- hT?«y. ^06lC 
/TxTcfcaiC, 7-yx$SMQSgpi 5tt> $*\ A7f 
^SQIISP 1 4 #f|=j$ L/i: A7J 0 - p ¥Ex-£ h 1 

(0, p) ZAtltZ (Xf7 76 0 1) o M 

smm^nicmmsMLr^k (X777602), 
FFT-rato-s*3S7-yxs«iic«fcy hi (0, p) 
icjj^Tpli73rp]<Di *7u7-yxs«i^t\ ^co/t 



(7) 
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7-*hi (0, p) iztmtZo osyjiift»p^p 



$SjH«H1 (0,p) £TIB (8) iCOCt<*46« 0 
[0 0 7 0] 



(7) 



(8) 



(Xt7^6 0 3) 0 

[0071] -f-LTs ns3me&>r>9v*>bi,rc 

#£5tK L Uf776 0 5) , 0 max*387?«S^ 
lcliXx i y76 o 3Stf 6 0 4dM0SHIUiEU W 

[0 0 7 2] CCT?ttA7J»6ii«l (x, y) 

lc*fjST3H1 (0, p) £fMBf*»&fcO^TWE 
Lfctf. #^«JSHfi!R (x, y) fcfcfjSfSHO 

(0, p) lco^Tt»±§S (7) S<OCi:<l^l«llc^fig 

[0 0 7 3] ±IB (4) s (7) s (8) S<fey HO 
(0, p) <fcH1 (0, p) tOBBfigttTIB (9) 
«fc?fc:SS*u (4) S£®il*frS¥ffl*»i®K»tf 

[0 0 7 4] 

H 1 (e,p)= 8 H 0 ($-t t sp) (9) 

#tfC, 02<DAx-y72O4 lcifxTftixJS*ISifl$!PI<D 
[0 0 7 5] ttSj£«$Sjft!g£-f*3Sttx 
[0 0 7 6] f&fo-Sx ^•i'J*s^t^TA7J;I 

sm»\ (x, y) ft^mmmw&R u. y) ic*jl 

H1 (0, p) &#PS 7- U Inflate HO (0, p) IC 
*fLTplA75ftlcf§/K ffi*L/ct<Dlc&oTl^ (0!l 
A7J»6ii«l (x, y) ##B3»£iiftR 
(x, y) ic^LTffiALTl^i*^ plA75ftlcf§ 

[0 0 7 7] £<D$$, ffi* • s Zm&tcOT* 

cl <t <L ^£^o 

[0078] z. x\ mmtipm*, m muDKWum 
aqaictt»j£«$«r*c£ic<fcy. pa^iRUcHfii 

#fff§ LTl^*IM», q*ll*lRllcli«^¥^j»Bi LT 
^SBBflHcSSjfSo Oty> #BB0-q¥®fc*fLT 
At: 0 - q ¥ffi£> q ii75[p]lC ffi* • f|g/J\$ s ICJS Cft 

[0 0 7 9] £0<fc ? IC q«l7J[RllC¥^» LTL^H 



[0080] m7 ^ Ltdmt&wmmmpym 

ixSISJUSQIIBP 19(i, $ 1 0 =&-tf □ fc«BKR£f * 
Ut7 7"7 0 2) 0 

[0 0 8 1] qfc*USf5p£TIE (2 5) 

p = c • exp ( q ) (25) 

q ±®*IIffiBlE*f JSf % p |i±(DJS*IfSB?-5R46^ 
Ut7 7"7 0 3) 0 

[0082] £5LTp<tq <D*f JSUfWUBJi Lfc&6 
«\ TIB (2 6) 5eicifx-fct:5^ 7-'JiSiI§H 

1 (0, p) »j»-r*7-uiM«jsaBE«iii«Hi 
(0, q) ics^r^o 

[0 0 8 3] H1 (0, q) =H1 (0, p) (2 6) 
Ut7 7"7 0 4) 

[0 0 8 4] 0$»A TIB (2 7) ^ICTxI-cfcaiC, A 
7777 0 4TW#*tlfcH1 (0,q) W\-><?Jm 
ixW (q) **J|Lfet©*, §?/c&H1 (0,q) t? 

[0 0 8 5] 

W (q) =0. 5 (1+cos (/rq/qmax) ) 

H1 (0, q) =W (q) • H1 (0, q) (2 7) 

fc/cU qmaxliqCDftAfiTS^ (Ax i y77 0 

5) 0 

[0 0 8 6] -D^T\ q£f>*y*>hL (Ax-y7 
7 0 6) s q #* ±m q max*3§T?£ fltf (At 1 V 7* 7 
0 7CDIUKYE S) , SffL/cqlc»LT^«WI 

(Ax-y77 0 3~Ax-y77 0 6) £1*yjlfjb\ •¥> 
^Tq^ftAfiqmaxlCjir^i (At V 77 0 7<DfU 
BtN 0) , 0^(DAx-y77 0 8lc^T?tl^o 

[0 0 8 7] ^-jsii, 06V>7y^>h?tis e#* 

*fi!0max*Si7»*H# (Ax-y 77 0 8(DIUBtY E 
S) % SiTL/'c0lC>PtLTq^BJS-lfPlcLfc±T« 
ftjflS (Ax'>77 0 2~Ax-y77 0 7) &m*m? 
t\ Wce&Wi*im6maxlcMT%£ (A7 1 'y77 0 



(8) 
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[0 0 8 8] ZZTMXjJimmm I (x, y) 

fcftfSfS H1 (0, q) £fMBf*»&fc:o^TKE 
Lfc*>\ #P,S;lJtl]#R (x, y) fcfcfj£f*HO 
(0, q) ICO^Tt»Mt5:t^'t5 0 
[0 0 8 9] C5Lfc^i*tlfcH1 (0, q) <tH0 
(0, q) OBBffitt, TIB (10) iCO<fe5[cS**l 

[0 0 9 0] 

tfi(0,<r) = 5.ffo{5-^ ? + A) (10) 

fc«U 

<Z = logp (11) 
A = logs (12) 

[0 0 9 1] C®<fc5fcA73»6ii« I (x, y) left 
J£f *H1 (0, q) li, #P,S;lJtl]#R (x, y) lC 
ftJfcfSHO (0, q) qli^r^lC-^, 01*75^ 

[0 0 9 2] ioL ^(Dqli*^(D->7hM-^^» 
ar«C<!:6^#ti^ ffi*-*i/h*s£±K (1 2) 
ie®i!Bffi<fcy$4&*c 0ii75ft<D~>7 

[0 0 9 3] C5L/"cqill*lRl05/7h«- 
A, 0|i75ft<D~>7 hmi/^&^C ^K)M»[iJgP 
1 6li, S*\ IWHx-^EfliaJI 7lc^46IB'ftLT33 

J -co Jo 



wdtmmmmimomffiT-* hq (0, q ) 

5) , co^*aj*ti/-c7-uxj5f8tjsa$Siii«Ho 

(0, q) tXf772 0 4*fMMrtlfc7-yxH* 
jS«$«ili«H1 (0, q) <t£ffl^T, Ztl64)2& 
TcfflBBfliaCr (0, q) U7772 0 

6) 0 ^<Df£, C®2 3»j7DiBBfliaCr (0, q) 
±<!:&3 0, qCDfeBcfcUHf^, ffi* s=& 
3<#>3 Ut772 0 7) 0 

[0 0 9 4] 0 8tt, C5L/-c2 3»j7ciiBBfliaCr 
(0, q) ®£JM!!g33<fclflH]£fc^ «**«'.h*s 
O»aj5BS0#|li€-^f 7P-f 1 h?* U <, IW1E18 
Kjfxfcfcdte, $rm (0, q) 52:^7-'Jlt 
JHLTF1 (u,v) £#463 Uf7 78 0 1) 0 oi/» 
T\ £<DF1 (u,v) ®&fiB#®/t7-£1. OKIEfc 
fbLT F1<p (u,v) ^46^ Uf7 78 0 2) 0 # 
PSJIJIH^R (x, y) KftjSfSHO (0, q) ICO 
^Tt^«cDSQII^HfT?tls COlSSfMB^flfc FO 
9 (u,v) ^IB'lSPlc^46IB'l?tlTl^o KT, 
£<DF0<p (u,v) #SI*aj**lS Ut778 0 
3) 0 

[0 0 9 5] ^^-»(DSQIIlco^TSlcf¥L<I/^ 
T3<t, $"f, HO (0, q) , H1 (0, q) 
^2 3*j7c7-yxS«S*tlTTE (13), (14) S 
lcifxTct:5lcF0 (u,v) , F1 (u,v) #£#6*1 

[0 0 9 6] 



Z.Z.X\ ±IB (1 0) iCtCtiS (1 3) , (1 4) iC 

fo ( U ,v) tFi ( U ,v) tcomim 

(15) =&||3o 
[0 0 9 7] 

= s ■ e-M*~W . F 0 (u,v) (15) 

FO (u,v) W7-ttiffl£Mf%t. TIB 
(1 6) ie®<fc5fcS**U F1 (u,v) fcO^Ttt± 
IB (1 5) iCcfcUTIB (1 7) iS<D<fcote^;*-*i£o 
[0 0 9 8] 

F 0 (ti,w) = |F 0 (u,i/)|c-J»*»(«.*) 



Fjfov) = Js J Fo(u,f)|e-J^- i ^ e -J^Cu,t ) 



(16) 



(17) 



d0dq (14) 



[0 0 9 9] 

foK») = e -**>(*.«0 



(18) 



(19) 



fCTs C5L/-cKfi**l/-cF0<p (u,v) , F1<p 
(u,v) =&ffll^TF1<p (u,v) £\ F0<p (u,v) 
*®j»£TK (2 0), (2 1), (2 2) j&Cifxf <fc 
5lCjfi7-UxS«iLTH0 (0, q) <tH1 (0, q) 
CD^lfi (2^*afll{£gO Cr (0, q) **d&* 0 
[0 10 0] 



£<DfcSJ|, FO (u,v) , F1 (u,v) 
tlftlTE (18), (19) 3i©<fc7fc*;!r*l*o 



(9) 
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oo J — oo 

<M**rt&xdo (20) 

= 6(0 -,(»,q -(-A)) (22) 

fcfcU F0<p (u,v) F0<p (u,v) 

Ut7^8 0 4) 0 
[0 10 1] C®<fc5fcLT$ttSftfciil!BfiICr 
(0, q) ttxVl/*l!B»fc&oT^*®#fofr* (±IE 
(2 2) 0 -f-LTs COxVl/^BBRCr (0, 

q) ^ 0 = <^ q = -A0>&flHft1T?ft*ttfr&3o 
Cr (0, q) ¥Ejb % 6«*fi!£<!:3SIf»B 
0 = 0, q = -A^ttJL (XT778 0 5) , £<D*t 
ffiffiB -A) lc*-5#< @Kft^^4&«ttt 
&**a/h*s«. _blB (12) iciest <$«i 

s = k -exp (2 8) 

6) 0 

[0 10 2] J-X±, C4>H8?f& HO (0, q) Icftf 
3H1 (0,q) <D S emiofa(D~S7 hw, qii£ft<D~> 
7 h«£#4&*/-c46fc\ 23»j7c7-UX$Si(3!>ffiffl«piJ 
ffl Lfc7- U Xffiffl$SiS6*ffl^S <fe a IC L/c6\ M% 
h7^/l/*^2:3^|!Bfc<fcoT*4&*||B6fc 

[0 10 3] :Mc, 02<DXx i y72 o sicifxf p*B5 
fflBBB«Cp (0, p) ®fMB#|iifco^T£ttWfcW 

[0 10 4] 0 9 ^ 0 2<DXx i y72 0 8®pfflSi 
BBHfcCp (0, p) (Dftmm\mm?7n-?^-h 

C p {9,Ap) = 

1^(6, p + Ap) 



[0105] mm9izmt&ofc. ®mm»&B-\ e 

it, m 8 lOTxTSQIUCcfc y ls]££ ffi* • |g/J\* s £ 
£4&fcfc£t*\ A^0-p¥Ex-£h1 (0, p) S 
tf^e-pTEx— £h0 (0, p) ^A^lLTs 3< 

totcEWS <i> JjcfctfffiA • iB/j\* s ^m^xmrn e-p 

¥Stc»r^A^l 0 - p¥E<D0li£ft<D->7 h«fi<fe 
tfpiA£ft<Dffi* • f§/J««IEir3 Ut779 0 

1) 0 

[0106] wz\ &me&e&-ecucmmBML 

Uf7 79 0 2) , reLMAplC-pmax X25ft 
AT3 (At'^9 0 3) 0 

[0 10 7] fLT, ffi* sic £o 

TSlEL/chO (0, p) <t hi (0, p) OIEfcfbiS 
fflBB«a*W-3tLTpiiSiBBli«C (0, Ap) lc|B 

«l U7779 0 4) s Ap*^r>^yy>h-rs 

(At'^9 0 5) 0 

[0 10 8] ^O'&s cKDAp^pmax X2*38T*S 
A^fr^StKL Uf779 0 6) , pmax x2*t^T- 
Xx779 0 4lc^TLT±IBAx-y7 
9 0 4Stf 9 0 SOiBS^Hyig-To 

[0 10 9] CtlfcfcfLTs ApftVmax x2J-X±<t& 
o/c^-g-lcli, ftJS$»0*>r>^yy>hL/-c^c 

ii-&lc«Axy79 0 3(C^-r« U7779 0 

8) 0 

[0 110] f&fo-S, ^O^K)M»[iJgP1 h 
0 (0, p) <th1 (0, p) (D#0lco^Ts pHi^IrI 

p^aiB^Cp (0, p) *f&*LTl**o 
[0 1 1 1 ] t6 Lm*M pT?S«»&0 1 *7U<DiE3i 

[0 112] 



-/ti)(Vfl,p)-/io)dp 



(23) 



v (x, y) <mmm^^xwm^o 

[0 1 1 3] 01 014, 02<DXx i y72 0 90ffi/\7 
£&SflH nv (x, y) ®fMB#|l!£^T 7P-^+ 

[01 14] wmi o^f<fcdic> ^KM»[ijgpi 6 

14. p*B5*BfHBiiCp (0, p) ZftmLtcKZlt 
U7771 0 0 1), fc«]JBR£L 

fcfcfc U7771 00 2) s plc-pmax 
(77771 0 0 3) 0 

[0 115] -f-LTs £<DCp (0, p) #RCp 
(0, p) ©ftAffiT&SCmax Jcyfc**!/^*^ 
(X77 71 004) s Cmax «fcy 



ICli, £<DCp (0, p) =g-Cmax leftALTCmax £ 

®rr z.<mo>pikp kicf-tAr^ (at 

7^1 0 0 5) 0 LfcAbT, £<Dp klcli, Cmax ft 

[0 116] :^lc, COp^V^yy^hL U77 
71 0 0 6) > p*Vmax *Sr?£*jb^fr£fltKL 

U7771 0 0 7) % pmax *3lT?S«»&lCl4. A 
7771 0 0 4~1 0 0 605BS^ijiyjIt"o 

[0 117] ^tllcJtfLTs pj^pmax lX±t^tcM 
Cp (0, pk) *ffi/\7$«i-r«o 0$»A 
5S!/\7$SiftlSSP1 8T?14. Cp (0, p) ±OjS 

(0, p k) £ x 



(10) 
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Pk — x<x&9 + ysm9 



[0 118] -f-LTs Z0MJ\7mtSlGk<D I n v (x, 
y) ¥E±<DlIII 

y=- (1/tan8) x + pk/s i n6 

±icc P (e, pk) (Dfi^o»r^ U777100 

8) 0 

[0119] -f-LTs nmm.e*'f>fij*>t>Ltc 
mz Ut'^i 009) „ modern* &mr*&%fr 

"Sfr&mmL Ut7 71 0 1 0) s 0max *tSTS^ 
ii-g-lcliXx-y^l 0 0 3lcf^Tf 3o 
[0 12 0] r^fe-Ss £<D^ftmiitiJBP1 6T1*s C 

p (e, p) ^e^eT-ft^fi^i^fSB (e, pk) 

llnv (x, y) ^J^LTl^o 

[0 12 1] ^LTs Z(Dm/\7%mW®. I n v (x, 
y) (DftAfiKDfeB (Xmax , Ymax ) #ff^ii 

(x A, y A) tK% 0 

[0 12 2] *3&Efc«fo*ii«BB^BB 1 0£ 

[0 1 2 3] 0 1 Hi, 0 1 lOTxI-HMS-g-gB 1 0£ 

[0 12 4] 011 (a) tt, X? ( r^j ) 

®#fiS»Sli«R (x, y) T?«y. 
R (x, y) 6^6x«y>>"*lRl^*4aLfc#BSX'y>>"*lRl 
IfRe (x, y) tt, 01 1 (b) lOTxI-cfcaic^: 

-So 

[0 12 5] ^-LTs c(DR e (x, y) KfcfLTM? 

$si*fifi-rt. 011 (e) izfrst&or^wcomm* 

JfO^e-pTSx— SihO (0, p) ££»A ?6 
lC£(DhO (0, p) £7-Ui$gjf 01 1 

(g) KaVf.fc^J&flfi&e-pTST-^HO (0, 
p) <Wi5*l«o ?6lC^(DH0 (0, p) Opft^q 
Hilc*f8tJS«$Si-r«t. 011 (i) lcifxTcfc5&# 
P,S0-q¥Ex-*HO (0, q) <WfS*l*o 

[0 126] 11 1 (c) tt, X? ( r^j ) (7) 

A^;1JSH# I (x, y) T-fe^o 

[0127] li^ft^tt^r^/cfeic, hicia 

z>ts £<D0i 1 (c) loTxTcfcaic, mfrtzntc*. 
mmmwmR (x, y ) ±<dx^icmlt> 0 

[0 12 8] C(DA7:;iJSij# I (x, y) fr6l-yv>" 
75^tttiJL/cA7JX-y5>^[R]iii# I e (x, y) tt, 
011 (d) IC/Txf £5K&3 0 

[0 12 9] ^LT\ CCD I e (x, y) lC*JLTM7 



011 (f) tcifxTct:5&A7J0-p¥ 
Ex-* hi (0, p) ££»A ?eicC0h1 (0, 
p) *7-'JJLmmt%ts 011 (h) fcifxf <fc-5& 
A7J0-p¥Sx-^H1 (0, p) jjftfS*l* 0 £6 
ICCCDH1 (0, p) (DpIi^qlilc^ixSIS^r^ 
01 1 (j) lc^-r<fe5SA*0-q¥®7 ? — 5»H 
1 (0, q) <WfS*l*o 

[0 13 0] £<Dh1 (0, p) #h0 (0, 

p) lcibLT#®5fay#jl,S*l*3Sttx 
ft<tffiB6^*Sfc46T*So 

[0 13 1] $/c, A7J;1JSB# I (x, y) 
J6H«R (x, y) left LTffiALTl^ /c#>lc, A7J 
7-Ul$gjij«H1 (©. P) (i#i7-'JlfiI 
fcHO (0, p) (C«LTpill7j(p)lC|iA/TVS03b^ 

[0132] $/c, zoLtcpibijmzwmimfoT^ 

ZMmts *f»JS«$Si**l/-cH1 (0, q) s HO 
(0, q) q»£lRKCH»«Wfc»LTV 

(0, q) ^T75lC->7h?-y:/ct(D6WH1 (0, 

q) ?65<:i6Wo 

[0 13 3] JT, 01 1 (k) tt, ffiUfiCr (0, 
q) ¥ST-Sys £<D*glfilCr (0, q) ¥®tts ffi 
OjStl4W5AScESiJT?#«**WJSOJS«ffiB*^L 

**HffiB#s 0 = 

q = - AT:&Z>) 0 
[0 13 4] $/c, 01 1 ( I ) 14. hi (0, p) £ 
hO (0, p) <D=&eiCOlv£\ ptt^lRllcrSLS^ 
5W-3tL/-cpfflSfflBBli«Cp (0, p) ^/TxT0T-S 

y, 01 1 (m) c P (0, P ) &e>&ev*±m 

^JfOffiB (0, pk) £fcajU ^cDfSBT-a4/\7^ 
«*fr-pfcjl6/\7$Sili« I nv (x, y) 
[0 13 5] ±5ELT*/ccfcaic, *HS6(DmJ8Ttts 
0 - q ¥®fMBSB 1 2 #I>y ^ffigp 1 3 s /\7$SISQ 

ssp 1 4 s 7- y xs^snusp 1 5 e<fetf!W8tjs«$Si 

iQSSPI 9*ffl^TA73B«lc*fJS-r«h1 (0, p) 
StfHI (0, p) £fMBU *SfcCft£*f»J£a@£ 
SILfcHI (0, q) «ft|*U ^ftmittiJBP 1 
<DH1 (0, q) tflB&T-^KttSn 7lc|B'lLfcH0 
(0, q) ^tbi$LT0-q¥®U^;l/T-0lE^^33cfc 
tfffi* • ffi/J\* s =&?R46^. <k <k & ic, fflsffiS H> fi<fetf 
ffi* • ffi/J\* s letter 0 -p¥iP^/R-¥W» 
M (x A, y A) *3ttB-r«<fe5«fiEL/-cOT?. 

ternr** =e y -y-<x^»$ij loo, b^cdx -y s^r^ji 

[0 13 6] 

[^<D^m] w±K^urc«fc5ic*«^-?tt, i\7m 

tA&fivIt 0 - p ¥S^* 5 ic 7- y i$» L/i: 0 - q 



(11) 
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0K£x ffi* • ffi/>*0f^£l»£Lfc±T?¥ffi&»* 
[0 13 7] -)Kfk/\7$SiA^5x 015^ 

[0 13 8] t'Jfeltv 0Kftd3<fclfffi* 

Rjsassir s c <t ic* y . #ps e - q ¥Eic*j lta 

73 0 - q ¥E6\ q IA75 ftlcffi* • ||/h* s U/cM 
[0 13 9] $tc, *£E®BiJ®38ET?tt. mmMz 

*^sp<fe y A7j-r«A7j#at*^^ #fisa«icBB-r 

A7:Bf»b^frSA7:T3fulc^4&§E'fI 
S*T£<«fc3fcLfc©Tx A7JB«^A7J**i/-cS5 
K#fiSB«t A7JB«t*IPJSfcfiS-&-r « C t^Rltlt 
& y , SQIIBiH^SlcJslir * C <k#T?#* <fc a left 

[HBOtBmSW^] 

[01] 01 te*mmmzwmz-£i&w(DmiM(DBm 

[0 2 ] 0 2 tt*&Efcft*B«SB^£5£®Hfifi®J&H 



[0 3] 0 3 1*0 1 icifxTi y 5*WHSMWt3x 5>* 
[R]B#0^J5K#|lI : &/Txr7P-?1'- hT'$5„ 

[04] 04ttv-^/i/©fs»^u-^^-r0if« 

[0 5] 0 5 1*0 1 fcifx-f /\7$«ft!gSftM55 0 - p 

[0 6] 0 6 »0 1 \zm 7- y x$si«isgi5^5 0 

[07] 0 7 »0 1 ic^-r*f»js«$Siftissp*^5» 
[08] 081*0 nc^-r»B»iwifflfl»fT3iiiisfta5 
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